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ABSTRACT 


A method of designing an optimal PI controller is discussed. 
The stability robustness and optimality robustness of the system 
using the optimal PI controller are discussed. The stability 
robustness is expressed in terms of the gain margin and phase 
margin, and in terms of bounds on system uncertainties which will 
preserve the stability of the perturbed closed loop system. The 
optimality robustness is expressed in terms of bounds on the 
perturbation in the system matrices such that the optimality of 
the closed loop system is preserved. 


Some numerical examples are solved to Illustrate the lesults. 
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CHAPTER - 1 


INTRODUCTION 

In this thesis the problem of designing a robust 

Proportional + Integral (PI) controller is considered and its 
stability robustness is Investigated. 

The servomechanism problem for systems with constant 
disturbances is a particular case of a more general problem which 

is considered by Davison [5] In his paper Davison has considered 

a more general class of disturbant “S and input for both regulator 
and tracking problem. The necessary and sufficient conditions 

were established for the existence of a robust controller for a 
linear time invariant, multivariable systems so that asymptotic 
tracking/regulation occurs independent of input disturbances and 
arbitrary perturbations in the plant parameters of the system It 
was shown that any such robust controller would consist of two 
devices (1) a sei vocompensator and (2) a stabilizing compensator 
The structure of the general robust ser vocompensator as suggested 
by Davison is shown in Fig. 1.1. 

The ser vocompensator is a feed back compensator with error 
input, whose dynamics depends on the disturbances and reference 
inputs to the system. For step input and constant disturbances, 
the ser vocompensator reduces to a integral controller. 
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The perpose of the stabilizing compensator is to stabilize 
the resultant system obtained by applying the servocomp ensator to 
the plant. 

The servomechanism problem has also been considered by 
Francis *£ clI . [8] by using a geometric approach. Some other 

papers dealing with the special cases of the problem discussed by 
Davison [5] are [6], [4], [18] and [20]. 

The stabilizing compensator stabilizes the closed loop 
system after applying the servocompensator . Such a compensator is 
not unique and may be designed in a number of ways e.g. by pole 
placement using state feedback [6] or by using ouput feedback 
[ 16 ], [ 21 ]- 

In [6] state feedback appt oach has been used and and 

(in Fig. 1-1) are calculated so as to place the closed loop poles 
at the desired location to achieve certain performance 
specifications. In [18] and [21] the problem of designing dynamic 
compensator, using output feedback for pole placement at 
prespecified location, is considered. This can also be used for 
designing stabilizing compensator. 

In [7] a compensator scheme, known as complementary 
controller, is suggested which can be used for the stabilizing 
compensator. In [24] output feedback Is used for stabilization. 
The determination of the required controller gain Is converted 
into the solution of an optimal control problem with quadratic 
performance index. 
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In this thesis state feedback Is used for stabl llzation . 
The problem of calculating gain K, and (In Fig. 1.1) Is posed 
as the linear quadratic regulator (LQR) problem. Here a quadratic 
cost function of the derivative of system states and the error 
(which In the case of a regulator problem is the output) and the 
derivative of control Inputs is minimised. The control function 
such obtained Is a linear combination of the states of the system 
and the Integral of the outputs of the system Is 

u X (t) = -K 1 x( t ) - K 2 J y(t ) dt 

The problem of designing the optimal PI controller can be solved 
as the standard linear quadratic regulator design problem. Then 
any of the algorithm available for the given regulator design can 
be used by replacing system matrix by augmented system matrix and 
Input matrix by augmented Input matrix 

In this thesis the designing of an optimal PI controller Is 
converted Into LQ regulator problem. Therefore the robustness 
properties of such PI controller can be related to the robustness 
properties of an LQR. Hence a survey of literature dealing with 
the robustness of LQ regulator is given here. 

Anderson and Moore [1] have shown that SISO liner quadratic 
regulators have infinite gain margin, 50% gain reduction tolerance 
and atleast 60° phase margin. Safonov and Athans [183 have 
studied the robustness of LQ regulators against large dynamical 
time varying and nonlinear perturbation In the f eedbackgains . 
They showed that If the control weighting matrix was diagonal the 
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LQ regulators would have infinite gain margin, 50% gain reduction 
tolerance and ±60° phase margin in each channel. This is an 
extension of the classical rsult due to Anderson and Moore £1 ] to 
iBUltilnput case. 

In [17] Patel et a.1 have expressed the robustness property 
of an LQR design in terms of bounds on the perturbations in the 
system matrices such that the optimal closed loop system remains 
stable. The bounds are expressed in terms of the weighting 
matrices. Katayama and Sesakl [12] have obtained several robust 
stability results derived for linear, nonlinear, time varying and 
dynamical perturbations in terms of bounds of perturbations and 
the weighting matrices assuming that the perturbations satisfy the 
matching condition. In [23] it is shown that the optimal state 
feedback law designed for nominal system can stabilize uncertain 
system provided the uncertainties satisfy the matching condition 
and are withxn given bounding set. Lehtomaki <et al [13] have 
derived robustness result of LQR using frequency domain criteria. 

OPTIMALITY ROBUSTMESS 

I t may happen that the opt lmal control law for the nominal 
system, in some cases, no longer remains optimal for perturbed 
system In view of some inevitable parameter variations in 
realistic systems, it might be more desirable and practical that 
the optimal control law remains optimal in presence of small 
parameter variations (optimality robustness). The obvious 
significance of such an optimal control law is that the associated 
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optimal regulator preserves the robustness properties in the 
presence of parameter variations. 

Hole [10] has shown that optimality is preserved under a 
given perturbation in the system matrices if the system is 
designed with a certain degree of stability. Barnett and Storey 
[3] have found the forms of perturbations in system matrix and 
input matrix which will leave the optimal (linear quadratic) 
control law unchanged. 

Fujii and Mlzushlma [9] have derived necessary and 
sufficient conditions for a multiinput optimal control with 
respect to a quadratic cost to remain optimal in the presence of 
small parameter variations. Mori and Kokame [15] have discussed 
necessary and sufficient conditions for the optimality robustness 
of a given control law when the system matrix is expressed as the 
interval matrix. In [23] also the robustness of the optimaillty 
of linear quadratic regulators is studied. In [22] it is shown 
that for a given optimal gain, the optimality of LQ regulators is 
preserved under small perturbations in system matrix and input 
matrix by modifying weighting matrices. 

The thesis is organised as follows : 

The Chapter 2 deals with the design method of optimal PI 
controller. In Chapter 3 the stability robustness of optimal PI 
controller in terms of gain margin and phase margin is discussed. 
Chapter 4 contains the robustness properties of optimal PI 
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controller in terms of the upper bound on the uncertainties which 
will leave the s/stem stable. Chapter 5 deals with the design of 
an optimal PI controller which has optimality robustness to the 
perturbations in system matrix and output matrix. Chapter 6 
concludes the thesis. 



CHAPTER - 2 


DESIGN OF THE PI CONTROLLER 


a. 1 I NTRODUCTI OM 

In this thesis the designing method of a robust PI 
controller is discussed. A number of methods are available for 
designing a PI controller. These methods include pole placement 
methods, converting the PI controller problem into an optimal 
output feedback problem or converting the PI controller problem 
into an LQR problem. In this thesis the last of the above 
mentioned methods, 1 . e . , converting the PI controller problem into 
an LQR problem is used because of excellent stability robustness 
properties of an LQR. 

This chapter Includes, the problem statement, the necessary 
and sufficient conditions for the existence of a solution to the 
problem of designing a PI controller, and the design method of the 
optimal PI cntroller. One numerical example Is also given to 
illustrate the design procedure 

2. 2 PROBLEM STATEMENT 

Ue consider a linear, time invariant system described by the 
equations 

x(t) = A x(t) + B u(t) + d i (21. a) 

= C x(t) + d 
^ x o 


y(t) 


(2.1b) 
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e(t) = y(t) - y r (t) (2.1c) 

vhere x(t) « R n is the state vector, 

01 

u(t) « R is the input vector, 

y(t) « is the output vector, 

y (t) « R p is the reference signal vector, 
r 

d , and d ate constant disturbance vectors, 

1 o ' 

and 

e(t) « R^ is the error signal vector. 

It is required to find a controller fot the system (2.1) 
such that the resulting controlled system is stable and the error 

e(t) ► 0 as t ► on for all x(t^) R n and for constant di and 

do and constant reference signal. 

Such a controller will be a Proportional + Integral (PI) 
type controller which is a particular case of the generalised 
controller structut e as suggested by Davison [5] 

2.3 NECESSARY AND SUFFICIENT CONDITIONS 

The theorem given below is applicable when the disturbance 
and the reference signal are constant. It is a special case of 
the generalised theorem given by Davison [5] which deals with 
the existence of the controller which will cause zero steady state 
error in the presence of a wider class of disturbances and input 
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reference signals. The modified version of the theorem as suited 
to our problem is stated here without proof. 

Theorem 2. i 


A controller which solves the problem stated above can 
always be found if and only if the following conditions are 
satisfied : 


(i) (A,B) is a stabllizable pair 

(il) (C,A) is a detectable pair. 

(ill) The number of inputs is greater than or equal to the number 
of outputs i e . m ^ p 


(iv) 


rank 


A B ‘ 

C 0 


n + p 


i.e. there is no transmission zero at the origin. 

Condition (lv) means that thet e should be no zero at origin 
for single input single output systems, otherwise the pole created 
at the origin due to the Integral controller will get cancelled by 
the zero at the origin. This will result in nonzero steady state 
error in the presence of constant disturbance. 

2. 4 DESIGN 

It is required to design a controller for the system (2.1) 
such that the resulting closed loop system is stable and the 

error e(t) * 0 as t » to Also the controller is to be designed 

by converting the original problem into LQR design problem. In 
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the LQR design the system states are penalised in the cost 
function. Hence in the following steps the error e(t) is 
incorporated into the systems states 


Rewriting equation (2.1) 

x(t) = A x(t) + B u(t) + d A (21. a) 

y ( t ) = C x(t) + d Q (2.1b) 

e ( t ) = y(t ) - y r (t) (2.1c) 

-> e(t) = C x(t) + d Q - y r (t) 

Let n.(t) = e(t) 

=> n 1 (t) = C x(t) + d o - Y r (t) (2 2) 

Combining (2.1) and (2.2) we get 


x(t ) 

II 

A 

O 

I 


'x(t) 

n 1 (t) 


C 

oj 


_n 1 (t)_ 


B ' 


■ d. ' 


o ‘ 


u( t ) + 

1 

- 


0 


d 

L o J 


- y r - 


or 


x 


A x(t) + B u(t) + d - 



(2.3) 


where 


A 


A 

C 


O' 

0 


~K is the augmented system matrix of dimension (n+p) x (n+p) 
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¥ is the augmented input matrix of dimension (n+p) x m 



Wow taking the derivative of equation (2 3), we get 


♦ ♦ ♦ 


X = 

A x(t) 

+ 1 u ( t ) 

(2.4) 

Defining 

x(t) 

♦ 

= X(t) 


and 

S(t) 

= U(t) 

(2.5) 


from equation (2.4) we get 


x(t) = A x(t) + B u(t) 

Here 

i r ^ (t) 

x(t) = x(t) = 

L V*) 

Thus by taking the derivative 
description (2.6) where system 


( 2 . 6 ) 



of (2.3) we get a new system 
state contains e(t). Also the 


constant disturbance terms vanished. 
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Now we will apply LQR design method to design a controller 
for system (2.6) which will make the system (2.6) stable in the 

sense that x » 0 as t * as. Such a controller will be a linear 

function of state x. 

Now v« take as the performance index (same as in the case of 
the LQR) 

a s 

J = J (x T Q x + R uj dt (2.7) 

o 

whet e $ is symmetric and atleast positive s emidef ini t e matrix and 
R is a symmetric and positive definite matrix. 

The u which will minimise J subject to (2.6) will also make 
the augmented system (2.6) assymptotical ly stable under the 
following conditions. 

[a] Pair [A, B] is completely stabillzable 

[b] The pair [A, DJ is completely detectable where D = Q . 

If conditions (i) and (11) of Theot em 2 for the existence of 
the solution to the PI controller problem are satisfied, 1 . e . , If 
(C,A,B) is stabillzable and detectable, the pair [A, B ] will be 
stabillzable [24]. This will satisfy condition [a]. 

The optimum u la given by [1] 

u = - R _1 B T P o X (2.8a) 


K x 


(2.8b) 
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~ -1 — T — 

where K = R B P^ and P q is the unique symmetric positive 
definite solution of the following algebraic Rlccatl equation 


X T P + ¥ A-P B R 1 B T P + = 0 

o o o o 


(2.9) 


From equation (2.8) and using the facts 


u = 


and 


x = 


x(t) 
e(t ) 


Ue get 



u(t) 

r "| x ( 1: ) 

= " K 1 K 2 

L J L ®c*) J 

(2.10) 

=> 

u(t) 

- - Kj_ i(t) - K 2 e(t) 


=> 

U(t) 

= - Kj x(t) - K 2 J e(t) dt 



Thus 

we get a PI controller where 

Is proportional gain 


and is integrator gain. The design method can be summarised in 
the following steps. 

Step 1 : Check for the controllability and observability of 

[C,A,B]. It should atleast be stabilizable and 

detectable . 

St«p 2 i Form augmented system matrix A and augmented input 
matrix B. 


where 

A 


'A 0* 

C 0 
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Stop 3 i Make proper selection of The weighting matrices R and Q 
5 should be selected such that the pair [ X ] is 

observable . 

St«p A i Solve the algebraic Riccati equation (2.9) to get "P . 

Use equation (2.8) to calculate K and partition it into 
[ ] as In the equation (2.10). 

2. 5 EXAMPLE 

-1.268 -0.045 1.5 952 

1.002 -1.96 8.52 1240 

0 0 -10 0 
000 -100 

0 0 

0 0 

Input Matrix B = 

10 0 

0 100 



Output Matrix C = 



It can be seen that [A,BJ is controllable and [A,C] is observable. 
Forming augmented system matrix A and augmented input matrix B we 
get 
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-1.268 

-0.045 

1.5 

952 

0 

0 


1.002 

-1.96 

8.52 

1240 

0 

0 


0 

0 

-10 

0 

0 

0 

= 

0 

0 

0 

-100 

0 

0 


1 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 


0 0 

0 0 

10 0 

0 100 

0 0 

0 0 

Taking R = and $ = I ^ . 

— —1 / 2 

It can be seen that [ A, Q ] is observable. Using LQR 

design algorithm optimal gain K is obtained 



K 


-.592869 .45601 

.599 7985 


.61191 3.06x10 3 -.803 

3 06x1 0 _ 2 4.763 .5958 


596 

805 


Hence proportional gain 


K 


1 


'-.592869 

.599 


. 45601 
. 7985 


61191 

3.06x10 


3.06x10 

4.763 


and Integral gain 



17 


.596 ‘ 

. 8025_ 

2.6 CONCLUSION 

The design method of the optimal PI controller is discussed 
and the necessary and sufficient conditions for the existence of 
PI controller are studied The design is done by converting the 
PI controller problem Into the standard LQR problem which enables 
us to use the standard algorithms available for solving LQR 
problem for calculating the proportional gain and the integral 
gain K 2 . 


- 803 
.5958 



CHAPTER - 3 


GAIN MARGIN AND PHASE MARGIN OF SYSTEMS WITH PI CONTROLLER 

3.1 I NTRODUCTI ON 

Gain margin and phase margin are frequency domain stability 
robustness measures . Anderson and Moore [1] have shown that 
scalar input linear quadratic regulators have a gain margin of 1/2 
to co and a phase margin of atleast ± 60°. An approach similar to 
that used in [1] will be used in the next section to evaluate the 
gain margin and phase margin of scalar input PI controller. 
Safonov and Athens [19] have extended the result due to Anderson 
and Moore [1] to multi input case. Anderson and Moore [2] have 
used slightly different approach to extend the result to 
multiinput case. This approach will he used to find the gain 
margin and phase margin for the multivariable case of PI 
control 1 er 

Finally some examples are given to verify the results 

3.2 GAIN MARGIN AND PHASE MARGIN CALCULATION 

For evaluating the gain margin and phase margin of the PI 
controller we assume that the system is free from constant 
disturbances because the gain margin and phase margin are defined 
for nominal system. 

Uniting the system equations (2.1) free from disturbances 


gives 
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x = Ax + Bu (3.1a) 

y = Cx (3.1b) 


Rewriting the design steps as discussed in Chapter 2 the augmented 
system is obtained as 


x = X x + E u 


where 



(3.2) 


and 


Minimising the cost function 
co 

J = J ^Qx + uRuj dt 
o 

with respect to u and subject to (3.2) gives the optimal control 
law as 


u 

where K 


- K 


IV 

X 


R 


-1 


B 


P 

o 


is the unique positive definite solution of the ARE 
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X T ¥ + ¥ X + $ - ¥ B R ' 1 B t T = o. 

0 0*0 o 

Starting irom the ARE it is shown in [1] that (see Appendix 1) 

£l+R 1/2 K £-j<oI - aJ B R' 1/2 J Jl + R 1/2 IC^j^I - Xj ¥ R~ 1//2 J > I 

(3-3) 

Scalar Case 

Consider a single input single output (SISO) system 

» 

x = Ax + Lu 

y = cx 

For SISO system equation (3.3) reduces to 



and R is a scalar. 

R can be put equal to one without losing generalisation. 
From (3 4) we get 
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For interpreting the equation (3 5) 


replace A by 


b by 


and k by [ kj k 2 ]. 


I-A 0- 

[c oj , 

[:] 


Then from (3.5) we get 


= > 1 + 


[ k i k 2 ] 

h [ 


ji»>I - A 0 


( j<^I - A) 

(joI-A)" 1 


JClf' 

LI fl 

j^j L o j 


2: 1 


= > 1 + 


h “2] ['I 

i 


( j^I-A) 1 b 


c(joil-A) * b 


> 1 


( k i + tH ( j 


ji.il - A| b j > 1 


(3 6) 


k 2 C 

For the interpretation of the term ( + - s — ) see Fig. 3.1 It 

J<>> k 2 c -1 

is obvious from the figure that the term ( k. + -s ) (j<jI-A) b 

1 

represents the open loop transfer function of the system with PI 


control let . 
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To satisfy equation (3.6) the polar plot of the open loop 
transfer function of the closed loop system will never enter 
the disc centered at -1+jO point and of radius 1. Since the closed 
loop system is stable (ensured by the design procedure), the 
Nyquist plot of the system will encircle the unit disc centered at 
-1+jO point u times m anticlockwise direction where u is the 
number of unstable poles of open loop transfer function. This 
ensures ± 60° phase margin and 1/2 to ® gain margin. This result 
is the same as that obtained by Anderson and Moore [1] for linear 
quadratic state feedback. 

Multivariable System 

For evaluating the gain margin and phase margin of 
multivariable systems the results given in [2] are used here. 


if 


In Fig. 3.2 if the closed system is stable for U(j<^) = I and 


[ 1 + V *H][ 1 


+ v 


H] * 1 


for all < a . 


([.] indicates complex conjugate transpose of the matrix [.]) 
then the closed loop system shown in Fig 3.2 will be stable when 


U* (j6i) + U(joo) > I 


for all <*i (for the proof see Appendix 2) 
Rewiiting equation (3.4) 


fl+R (-J.I - s]'Ss lf2 l(j«I - a]' 1 b r- 1 ' 2 ] > I 
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which can be written as 

J^I+^R 1/2 K £j<oI - aJ I R _1/2 | J £l+R 1/2 K^jwl - aJ B R~ 1/2 J> I 

Replacing A by 

b by 

and K toy [ Kj K 2 J, 
we get 

[* + {* 1/2 [ K i * ] H - A ] _1 B }‘] 

[l * {r 1/2 (K, * ^ ) (j-I * A) \)] > I (3.7) 

Redrawing the Fig 3 2, taking 

V - R *'* («, ♦ V J {juI . 

and 

U = R 1/2 L R _1/2 

Ue get Fig. 3.3. 

Using the result discussed earlier, the system shown in Fig- 
3.3 will be stable when 

R' 1/2 L* R 1/2 + R 1/2 L R -1 ^ 2 > I (3.8) 





Fig 3 3 



Fig 3 4 System with PI Controller with L(jw) representmg perturbations 






27 

\ 12 . 

Multiplying above equation by R ' from both side. Ue get the 
equivalent condition 

L* R + R L > R (3.9) 

Nov suppose R is diagonal and L is drag (i.». ,1 ) then equation 

i in 

(3.9) is satisfied if and only if 

l* + > 1 ( i = 1 , . . , m) (3.10) 

If i. is real, then (3 10) will be satisfied for i < l < co. This 

i 2 l 

implies that there is a gain margin of i to oo m each channel. 

Also, if = e ^ , then (3.10) will be satisfied for ^ . 

This means that there is a phase margin of ± 60° m each channel. 

It can easily be seen that the two systems shown In Fig 3.3 
and in Fig 3.4 have same open loop transfer function matrices. 
Hence the two systems have same stability properties. Therefore, 
the results derived above for the system shown in Fig 3.3 are 
also valid £oi the system shown in Fig. 3.4. Hence the system 
shown in Fig. 3.4, also have a gain margin of i to <*> in each 
channel and a phase margin of ± 60° m each channel, which is a 
multivariable system using the PI controller designed using LQR 
design method. 

3.3 DISCUSSION 

There are a number of design methods available for PI 
controller. In this thesis a design method is chosen in which the 
proportional and integral gains are calculated by converting the 
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PI controller design problem into an LQR design problem for an 
augmented system. Now it is shown that the system with PI 
controller designed in this way shows a phase margin of atleast 
60° and a gam margin of to oo, same as shown by an LQR. 

3 . A EXAMPLE 

Here 2 examples are solved to verify the results proved in 
this chapter. 

0 1 01 

0 0 1 

2 3 4j 

' 1 

Input vector b = 0 

0 

and output vector c = £ 1 1 1 



Taking y = I ^ and R = 1 

The optimal proportional gain = j\l.694 21.516 34.047 

and the optimal proportional gain = 3.1381 

Calculating the open loop transfer function (GH), we get 



11.7 s 3 + 18.16 s 2 + 14.21 s + 3.138 



2s 
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GH has one unstable pole. The Nyquist plot shown 
encircles the unit disc centered at -1+jO once in 
clockwise direction as expected 
2 . 


A 



1 0 

0 1 

3 4 


b = 


' 1 ' 
0 
0 


and 


C = 




1 


4 


Taking Q 


1 0 0 O' 

0 0 0 0 

0 0 0 0 

0 0 0 0 


and R = 1 , 


the optimal proportional gain, we get 


k 


1 


£ll. 9 


-21.3 


-34 


] 


and the integral gain 


k 


2 


in Fig. 3.5 
the counter 


-.6 



■n<G(ju)rUjv)) 



I m(G(jw)H(jw)) 



-10000 -8000 -6000 - 4.000 -2000 0.000 

Re (G(jw)H(jw 

Fig.3 6 Nyqmst Plot of (11 9*s^+19 8*s 2 + 




The open loop transfer function (GH) 


K C -1 

GH = [ K 1 + “i-j( sI - A J b 

™ - li-9 s 3 + 19.8 S 2 4 10.5 s 4 3 

— / UH — T X X 

s - 4s - 3s 4 2s 

GH has two unstable poles. The Nyquist plot (Fig. 3.6) of th. 
system predictably encircles the unit disctwice in the count 
clockwise direction. 

3. 5 CONCLUSION 

The stability robustness of the PI controller in terms 

gain margin and phase margin is studied. It is shown both t 

scalar input case and for multiinput case that the gain margin 

-i- to oo and the phase margin is of atleast ± 60°. Some examp! 

<6 

are also solved to verify the results. 



CHAPTER - 4- 


ROBUSTNESS OF THE OPTIMAL CONTROLLER IN THE PRESENCE OF SYSTEM 

UNCERTAINTY 


4. 1 INTRODUCTION 

In practice, a system is usually subject to parameter 
variations and modelling errors which are often very large and 
which do not allow for an accurate mathematical representation of 
the system. It is therefore necessary to investigate the 
stability robustness of the PI controller design in the presence 
of system uncertainty. 

The stability robustness of an LQR in the presence of system 
uncertainty has been studied by Patel et al . [17]. In this 
chapter the robustness property of the PI controller is expressed 
in terms of bounds on the perturbations (modelling errors or 
parameter variations) in the system matrices such that the closed 
loop system remains stable as done in [17] for an LQR design. 
Relations are established between allowable perturbations and the 
weighting matrices in the quadratic performance index, thereby 
helping to select appropriate weighting matrices in the quadratic 
performance index to attain a robust design. Also the stability 
robustness of the PI controller is compared with that of an LQR 
design by comparing the bounds on the perturbations for stability. 

Finally some numerical examples are given to illustrate the 
results derived in this chapter. 
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4.2 PROBLEM STATEMENT 


Consider a linear time invariant system whose dynamics are 
expressed by state equations given below 

x(t) = A x(t) + B u(t) + £ [ x( t ) , u(t)] (4-la) 

y(t ) = c x(t) + g[x(t)3 (4.1b) 

where x(t) <e R n , u(t) « R m , y(t) <e R**. 

A, B and C represents the system model parameters, and £ and g 
nonlinear vector functions. 


Let y^ be the constant reference signal vector then e(t) = 
y(t) - y^ as defined in Chapter 2. 

Using (4.1b), we get 

e(t) = C x ( t ) + g [x(t ) ] - y r 

Following the design steps discussed in Chapter 2 we get the 
augmented system equation 


X(t) 


*f [x(t), u(t ) ] 


A x(t) + B u(t) + 


£ [x(t)] 


where 


x(t) 


r X(t) ■ 

L e(t) J, 


(4.2) 


u(t) 


u(t ) , 
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A = 


and B = 



In Chapter 2, the LQR design method was used for the design of PI 
controller. Here a slightly modi fed cost function is used to 
design the PI controller. 


The performance index to be minimised is 

Ob 

J = ^(t) Qx(t)+u t (t)Ru(t)Jdt (4.3) 

o 


where Q and R are selected as in Chapter 2 and a is a real 
positive number and is used to prescribe the degree of stability 
of an LQR design. The optimization of the performance index (4.3) 
subject to the system model 


x ( t ) = A x ( t ) + B u(t) 


(4.4) 


yields the control 


u(t) = - R' 1 B* P^(t) (4.5) 

where F >0 is the solution of following algebraic Rlccatl 
equation 


(a ♦ «.] 


ot 


+ p 


OL 


£a + otlj - 


P B R 

ot 


-1 


B 1 


a 


+ Q = 


0 


(4.6) 
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The resulting closed loop system is given by 

' f [x(t)] ' 

♦ 

X(t) = (A - B R -1 I* P^Jxft) + (4.7) 

£ [x(t) ] 

and is known to be stable, if 
' £ tx(t)] 

. £ [X(t)] 

The problem to be investigated is to determine the bounds on 
f x(t) 

& X(t) 

which preserve the stability of (4.7). It is the stability of 
augmented system which ensures zero steady state error. 

For finding the bounds a theorem given by Patel et al . [17] 

Is used here. The theorem given in [17] is applicable for LQR 
design. The theorem is applicable here also because PI controller 
design problem is converted into LQR design problem. The proof of 
the theorem is given in the Appendix 3. 

Before stating the theorem certain definitions and notations 
are given here. 

Euclidean norm : || w ||g (||U|| E ) denotes the Euclodean norm of 

vector w (matrix U) . 

rv, 21 1/2 

II“He“ [2 l u iji ] 
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x X(U) denotes maximum eigen value of U and min X(U) denotes 
nimum eigen value of U. 

ectral norm : ||U Jj^ denotes the spectral norm of matrix U 

ll°H. * -ax [x (u« T )] 1/2 = „ ax [x [« T 0]] 1/Z . 

K>or*m A. 1 i Let D = Q + ¥ BR^'B t P. Then if, the 

ca a 

mllnear vector function 

r f x(t) 

[ « xct) 

atisfies the condition 


f(x) 

fi(x) 


min X(D) a min XCP^) 

2max X(P ) max X(T? ) 
or - or 


for all x <e R 


n+p 


(4.8) 


the closed loop system (4.7) is stable. 

-inear Per tur bat. ions 

Now a special case of the perturbations f [x(t), u(t)) and 
g [x(t)J will be considered. It is assumed that £ [x(t), u(t)] 

is linear in x(t) and u(t) and g [x(t)) is linear in x(t), and can 
be written as 


f [x(t), u(t)) 


E x ( t ) + F u(t) 


(4.9) 



38 


6 [x(t)] = G x(t) (4.10) 

where E, F and G are nxn, nxm and pxn matrices r especti vely. 

The system equations (4.1) now becomes 

x(t) = (A + E) x ( t ) + (B 4 F) u(t) (4.11a) 

y(t) = (C 4 G) x(t) (4.11b) 


The matrices E, F and G may represent modelling errors or 
parameter variations in the system model (A, B, C). 

The augmented system equation (as in Section 2.4) now 
becomes 


♦ ♦ _ 
X (t) 

fA+E 01 rx(t)l 1 

"B 4 F" 

♦ 

+ 

n 


e(t ) 

I.C+G OJ j 

0 


u(t ) 


(4.12) 


which can be written as 


'x (t) 
e(t) 


:] [::::] • [ 


"E 

G 


U(t) 4 


0] jx(t)' 
Oj [e(t) 


[ 3 * 


(t) 


(4.13) 


or 


X ( t ) = X X(t) 4 E u(t) 4 "F u(t) 4 1 x(t) 

where x(t), u(t), A and B are defined as in Chapter 2 and 


(4.14) 
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'E O' 

G 0 . 


lomparing (4.14) with (4.2), ue get 


f tx(t), u(t ) ] 


l & (x(t)] 


B E x ( t ) + F u ( t ) 


Using the control law (4.5), (4.15) becomes 


f tx(t), u(t)] 

£ [x(t ) ] 

L 

where 


( E+ F K^) x(t) 


(4.15) 


(4.16) 


K = - R _1 I* P (4.17) 

a a 

Using the same control law (4.5) in system (4.14) the closed loop 
system becomes 


x(t) = 


A + E + Tb + f] K 

L l J <*J 


x(t ) 


C4.18) 


The following theorem which is derived from Theorem 4.lglves bounds 
on E, F and G for the stability of closed loop system (4.18). 
Theorem 4,. 2 s If the perturbation matrices E, F and G satisfy the 


condition 
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mi.* imui 


ot 



min X(D) ot min X(P^) 

2 max X(P ) max X(P ) 
ot ot 

(4-19) 


then the closed loop system (4.18) is stable. 


Proof 

From (4.16) 


|f(X(t), u(t)|| 

IgCx(t)) He 

ll*(t)ll E 


J j CE + F K^) i(t)jj E 
I l E 


< max 
x(t) 


I I (E ♦ r K a 3 *(t)|| E 
l|5(t)|| E 


or 


or 


* n e + f ui. 

* imu + imi. ii*ji„ 



Now using Theorem 4.1we get stability condition 



K 

ot 


min X(D) ot min X(P qi ) 
2 max X(P ) max X(P ) 

Ot ot 


Thus the theorem is proved. 
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Theorem 4 . 3 s If the perturbation E and F satisfy the condition 


'E 

G 





* P 
s 


min X(D) <* min X(P ) 

+ 

2 max X(P ) max X(P ) 

HOC ot 


C4.20) 

then the closed loop system (4.18) is stable. 

Proof 


The result of Theorem 4.3 follows from Theorem 4.2 on noting 
that for any matrix U, | jU| | > I |Uj | and therefore 

IL & * 



4.3 COMPARISON OF ROBUSTNESS OF PI CONTROLLER WITH THAT OF LOR 
FOR LINEAR PERTURBATIONS 
Consider system represented by (4.11a) 


x(t) = A x(t) + B u(t) - 4 E x(t) + F u(t) 


Let a linear quadratic state feedback (LQSF) control law be 
designed for the nominal system 

x(t) = A x(t) + B u(t) 

The performance index to be minimised is 


oo 


‘ J [ 


x(t) Q x(t) + 


u t (t) 


R u(t) I dt 


(4.21) 


o 



42 


where Q £ 0 and R > 0 

Assuming that (A, B) is controllable and (A, Q 1 ^ 2 ) is observable 
the optimization of (4.21) yields the constant gain state feedback 
control 


u(t) = - R 1 B* P x(t) 


(4.22) 


= K x(t) 

at 

where K = - R -1 B* P . 

o i at 

Uhere P > 0 is the solution of the following ARE 

fA + all P +P Ta + all - P B R _1 B* P +0 = 0 (4.23) 

l. J a a ( Ja a 

The resulting closed loop system will be 


x(t) = 


r r i 

|A + E + B + F 

L J 


a 


x (t) 


(4.24) 


Patel et al [17] have found bounds on perturbation matrices E and 
F for the stability of the closed loop system (4.23). These 
bounds are given by the following theorem. 


Theorem 4. 4 : If the perturbation matrices E and F satisfy the 

condit ion 


imu + 



min X(D) a min A(P ) 

< ^ = + — 

2 max \(P ) max x(P ) 

ot at 


(4.25) 


then the closed loop system (4.23) is stable. Here 
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D = Q + P BR 1 B t P 

ot a 

For comparing the stability robustness of the PI controller with 
that of LQSF controller in the presence of linear perturbation it 
is assumed that the perturbation matrices E and F are same for 
both the cases. Also the input weighting matrix R is same for 
both the cases. The state weighting matrix Q is different from 
and they are related by 


0 = 


Q 

t 


t 

s 


where t and s are chosen so as to make (j atleast n.s.d 

For comparing the robustness of PI controller we will use 
results of Theorem 4.3 and Theorem 4.4. First of all the left 
hand sides of equation (4.20) and equation (4.25) are compared. 

From the definition of Euclidean norm 


E 

G 



(equality holds when G = 0) 


and 




(4.26) 


(4 27) 



Now it is required to compare IJK^JI with 
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I l*J I 8 “ [* a *«]] 

Putting K = - R -1 B* P 

Ot CA 

and using the fact that 


max 





( 4 . 28 ) 


Ue get 


IlSjIa = ffiax 


K r „ b R ' lR " lEt p j] 


1/2 


= O. [p a B iT 1 ] 


(where ct [.] represents maximum singular value of [.]) 

-> H*Jl. * s [ R J 5 [ E E_1 ] 


or 




M s {[ :] ■-} 


since B = 


-m 


or 


< a 




(4.29) 

(4.30) 


Hence 

ll*JU 4 5 [ R J 5 [ B R ' 1 ] 
as 5 [ B 0 ] = 5 [ B r1 ] . 


(4 31 ) 
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Similarly 


11**11. - «“ (*I **]] «- 32 > 

=> ll**ll. - 5 [*I] 

Putting K = - R _ 1 B t P we get 

a ot 

ii**ii. ■ * [ p * b r_i ] 

j> n**n. s 5 [ p *] 5 [ b * _i ] <«- 33 > 


It is obvious from equations (4.31) and (4.33) that to compare 

||K^JJ q and ||K^||, it is required to compare o- [P^] There is no 

direct relation between o* [ P ] and o- fP ] since they are the 

at at 

solution o£ two different ARE ’ s (4.6) and (4.23) respectively. In 

literature, bounds on the maximum eigen value (which is the same 

as maximum singular value as P and ¥ are symmetric and p.d ) of 

ot at 

the solution of ARE are available. These bounds give an 

indication about which of o* [P ] and o [P ] is greater 

ot at 0 

Using one such result given by Mori and Derese [14] we get 

/ Z Z Z — LT'-t" ' 

tr (A + otl) + V tr (A + otl ) + tr(Q) tr (B R B ) 

tr (B r” 1 b 1 ) 

(4.34) 



also 
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" [ p «] s <*1 = 


tr (A + al) 


+ oil ) + V tr (A + 


(A + oil) + tr (Q) tr (B R 1 B* ) 


tr (B R 1 B* ) 


Here tr [.] represents trace of matrix [.]. 


Comparing ft and ft term by term 


(4.35) 


tr |~A + atljj 


tr 


A+otl 0 


°J y tr i k + “ : ] 


(4.36) 


tr { Q j > tr Q 


(4.37) 


aa y is of 1 at get dimension than Q . 

(Here it is assumed that for comparison purpose Q is chosen 
such that it contains Q in it i . e . Q is larger form of Q ) . 


tr 


[b R ' 1 B*] = tr |[ ®] R " 1 [b* 0 ]} 


= tr 


{[ ■ :"] [■• •]} 


= tr 


ffB R 1 B* 0] 


it 


0J 


= > tr 


|b R -1 B*] = tr [b R 1 B^ 


(4.38) 


From inequalities (4.36), (4.37) and (4.38) it follows that 
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fit > ^ 


Also, £pJ £ ft 1 


and 


W a r, > 


(4.39) 


(4.34) 


(4.35) 


Though It la 

not possible to say any 

t h i ng 

with 

certainty 

about 

the relation 

between <y [P ] and 

[P 3 

a 

but 

the above 

three 

relations do 

give an indication that 

& [P ] 

tends 

to be , 

greater 


than cr [ P ] . Hence similar kind of statement can be made about 

a. 

ll^oJIs an< * | 1 1 i s from equations (4 31) and (4.33) that I 1^*1 I a 
tends to be greater than IjK^JJ . This is verified by numerical 
examples also. Using tills and equations (4.26) and (4.27) it can 
be stated that 



Thus the stability margin in the left hand side of equations 
(4.20) and (4.25) is smaller in case of PI controller as compared 
to that of LQR. This Implies that PI controller is inferior to 
LQR as far as stability robustness to the perturbation in system 
matrices is concerned. It is the price that is to be paid for 
constant disturbance rejection and zero steady state error which 
are not guaranteed by an LQR. 

Now we will compare right hand sides of equations (4.20) and 
(4.25), i e., p and p. 
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nun A. ( D) ot min X(P ) 
ot 

2max X(P ) max X("P ) 

ot ot 


(4.20) 


D = Q + P BR 1 B t P 

ot ot 


and 


min X(D) ot mm MP ) 

+ ±- 

2max X(P ) max x( P ) 

CA Ot 


(4.25) 


In the literature bounds are available on the min. singular (or 
eigen) value of the solution of ARE. But these bounds depend on Q 
and Q and even the bounds on min X (P ) and min x(P ) can not be 

ot ot 

compared but it is seen from numerical examples that min X(]?^) is 
smaller than min X(P^). Now snin ^(D) and min X(D) will be 
compared . 


min X(D) = o* ID as D is symmetric p.s.d. matrix 


or 


■ 


Q + P B R 

ot 


1 B J 


or 


« \Q\ + & 


j^P a B R 1 B* pJ 


or 


2 ?[«]*? [ P J £ [ B * _1 B '] £ [ P J 


or 


> a 


[5] * £ [ F J £ 


B R 1 B t 0‘ 


-H 


= > mm X(D) > 


£ [o 


(4 40) 
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Sl&Ilerly it can be seen that 


mm A.(D) > </ 

Q 

+ <y fp 

1 |"b R~ 1 B***! 

- kl 

C 4 41) 



- L 

“J - L J 

- L «j 


F roc the structure of Q and Q it is obvious that o* [Q] > <r[ Q ] . 

Though nothing can be said about min X ( D) and min X(D) with 

certainty but inequalities (4.40) and (4 41) gives an indication 
tha mm x (D) is larger than min x(D). In most of the cases the 
difference between max X(P ) and max X(P ) is so large that 

CA OL 

relation between ju and Jj is mostly determined by max X(P^) and max 

X(P ). And since max X(P ) tends to be greater than max X (P ) 
or ol or 

(from our earlier discussion), p tends to be smaller than p is 
most of the cases. Thus the stability margin is smaller for the 

PI conti oiler than for the LQR . This is the price to be paid for 

constant disturbance rejection. 


4.4 HUMERI CAL EXAMPLE 


The results of numerical examples are presented in the form 


of tables. A 

llnearizedd model o£ a 

gas 

turbine 

considered 

Uang 

and Hunro 

[24] is given 

by 






-1.268 

-0.045 

1 . 5 

952 


0 

0 


1.002 

-1.96 

8.52 

1240 


0 

0 

x = 

0 

0 

- 10 

0 

x + 

10 

0 


0 

0 

0 

-100 | 


0 

100 


_ 



J 



- 
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y = 


1 

0 


0 0 
1 0 


O' 

0 


x 


Table 4.1 gives the results for various values of Q and (J 
for the above system 

Table 4.2 gives the results for various values of Q and Q 
for the following system 



M 

O 

O 

I 


1 

1 

* 




X = 

-10 2 

x + 

0 


0 

1 

»-* 

I 


O 

L. 


y = [1 1 1] x 


ol - .2 and R = I 

£ 


LQR design 

PI controller design PI controller design 



i 4 o] 



For Q = I 4 

For Q = 

0 0 

For Q = 

1 6 

1 l K «l U " 3 ' 69 

I|KJ»S = 

4 .27 

nails' 5 

. 278 

W D > = 1 

Snin^ ‘ 

= .007 

Nnm^ = 

1 

X . (P ) = .0015 
mm <* 

X . (P ) 
mm or 

= .0012 

X . (P ) 
mm or 

= .0012 

X. (P ) = .226 

max or 

X (P ) 

max or 

= .735 

X (P ) 

max or 

= 2.228 


Table 4.1 
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d = . 

25 and 

R = 

1 


LQL design 

PI controller design 

PI controller design 

For Q = I 

3 

For Q = 

1 3 °' 

0 °. 


For 0 = 

*4 

limits = 

= 12.81 

■ 

= 17.43 


IIKJL 

= 25.93 

^n (D) = 

: 1 

X mn (5) = 

= 0 


X min^ 

= 1 

CP ) 

m j n ol 

= 2 43 

\un Cl V 

= 421 


X (P ) = 1 .193 

nun or 

K (P ) 

max ot 

= 100 23 X (P ) 

max or 

= 168.621 

X (P ) = 324.99 

max or 


'l 1 1 


"ill 

0 


"l 1 1 l" 

For 0 = 

111 

For Q = 

111 

0 

For $ - 

1111 


111 


111 

0 


1111 




0 0 0 

0 J 


1 1 1 1_ 

H K JU : 

= 15.09 

H*JU : 

= 20.55 


IlSjl. 

= 25.06 

3 

W D) : 

= 0 

W®> 

= 0 


W°> 

= 0 

X . (P ) 
min oi 

= 2.75 


= 719 


X . (P ) = .381 
mm or 

^ CP ) 

max oi 

= 136.48 X (P ) 

max or 

= 225 32 

X (P ) = 310.83 

max or 


Table 4. a 


3 140b,- 
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From the results it can be seen that all the examples 
discussed in this chapter suggests the following. 


K | | £ 

CA* 1 S 


[K 


Ol * 1 Q 


X (D) £ X . <D) 

mm mm 


x (P ) 5: X (P ) 

mm or nun or 


and 


a (P ) £ X (P ). 

max or max or 


The above inequal i t i es imply smaller stability margin for PI 
controller as compared to that of an LQR design. 


4-5 COMCLUSIOM 

In this chapter sufficient conditions have been obtained for 
the stability of the optimal PI controller in the presence of 
linear and nonlinear perturbations The robustness is expressed 
in terms of bounds on the per tut bat ions in the system matrices 
such that the closed loop system remains stable. 

A comparison has been made between the stability robustness 
of a PI controllei and that of an LQR design Some numerical 
examples are given to illustrate the results given in this 


chaptei . 



CHAPTER - 5 


OPTIMALITY ROBUSTNESS OF THE PI CONTROLLER 

S. 1 I NTRODUCTI ON 

It has been shown In Chapter 3 that the optimality conditon 
(3.3) is responsible for excellent robustness properties in terms 
of gain margin and phase margin. But these properties are for 
nominal systems and can not be guaranteed for perturbed system. 
It Is desirable to have these properties for perturbed system also 
for which the fulfilment of optimality condition is required for 
the perturbed system. 

Hole [10] has given a sufficient condition for an LQR design 
which when satisfied guarantees the satisf act luon of the 
optimality condition for the perturbed system using nominally 
optimal control law when fixed perturbations occur in plant state 
matrix. In this thesis similar result is given for the PI 
controller . 

Hole and Hahanta [11] have given a sufficient condition for 
an LQR design, the satifactlon of which gurantees specified gain 
margin and phase margin for the perturbed system, using the 
optimal control law obtained for nominal system, in the presence 
of perturbation in state matrix. Similar result also exists for 
the PI controller and is given In this thesis. 



A bound on the perturbations in the system state matrix and 
output matrix is also obtained which will preserve the optimality 
of the PI controller under perturbation in the system matrix and 
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in the output matrix The applications of the results are 

demonstrated through numerical examples 

5.2 PRELIMINARIES AND PROBLEM FORMULATION 

Consider a linear time invariant system described by the 
state equations 


x(t) = 

(A + 

AA) 

x(t ) + B u(t ) 

(5.1a) 

y(t) = 

(C + 

AC) 

x(t ) 

(5.1b) 


where x lR n , u e fR m and y <£ . 

A, B and C are the nominal state matrix, nominal input 
matrix and nominal output matrix respectively, and AA and AC 
represents perturbations in the system state matrix and in the 
output matrix respectively 

The corresponding nominal system is described by 

x(t) = A x(t) + B u ( t 1 (5.2a) 

y (t ) = C x(t ) (5.2b) 

Now to design the optimal PI controller the nominal 
augmented system equations obtained as discussed in Chapter 2 are 


given below. 
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x(t) = X x(t) + I 0 ( t ) 


(5.3) 


where x 



y r (t) (The reference signal vector) 


and 


u - u . 

Similarly perturbed augmented system is described by 


x = X x + E u 
P 


where A 

P 


A+AA O' 
C+AC 0 


B 



(S 4) 


Augmented perturbed system matrix A^ and augmented nominal system 
matrix A at e related by 


A = A + AA 
P 


where 


AA 


AA O' 
AC 0 


(5.5) 


(5.6) 


For designing the PI controller the following cost function 
is minimised with respect to u, subejct to (5.3) 


<n 


J 


2 ott f~t K ~ ~ t ~l 

e x Qx+u Ruj 


dt 


(5.7) 


o 



where ot Is real positive scalar, 
Q > 0 and R > 0. 
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Under the assumption that the pair [A, 1] is controllable 

— —1 / 2 

and the pair [A, Q ] observable, a unique optimal control is 
given by 

u = -K x = -R -1 B 1 P x (5 8) 

a ot ' 


P^ in (5.8) is obtained as a unique symmetric p.d solution of the 
algebratic Rlccatl equation given below 


(I ♦ «.] P a ♦ (X * «>] - P„ B p- 1 s* f, * 5 

Using (5.8) in nominal system (5.3), we get 

x(t) = - B R~ 1 B* Pj x(t) 

It is shown in Appendix 1 that the closed loop 
satisfy the optimality condition given below. 


= 0 ( 59 ) 


(5.10) 

system (5.10) 


£l+R 1/2 £>1 - aJ 1 B R 1/2 J [i+R 1 /2 K^jwl - Aj B R 1/2 j > I 

(5 11) 

In Chapter 3 it is shown that the satisfaction of (5.11) ensures a 
phase margin of atleast ±60° and a gain margin of to ® for 

nominal closed loop system (5.10) using optimal control law. 

Now the problem can be stated as follows. 

It is required to find a condition under which the 
optimality condition is satisfied by the perturbed system (5.1) 
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using the control law designed for the nominal system, i.e , to 
find a condition under which the perturbed system will satisfy the 
following condition. 



(5.12) 


The satisfaction of (5.12) will ensure a phase margin of 
atleast 60° and a gain margin of ~ to oo for the perturbed system 
using control law designed for nominal system. 

5.3 OPTIMALITY ROBUSTNESS RESULTS 

The results are given in the form of the following theorems. 


Theor em 5.1 


Assuming that the pair [X^, B] is controllable and the pair 
— —1 / 2 

[ Ap, Q ] is observable, the control law (5.8) will satisfy the 


condition (5.12), if 


A - B K ] P + P Ta - B K ] + 

P a) a p 


K t R K 
at at 


(5.13) 


is atleast n.s.d. 


Note : It can be seen that 


F = AA P + P AA - IQ + 2otP 
a a a I at 


(5.14) 


Proof 

The proof of the above theorem follows on the same lines as 
given by Hole [10] for a similar result for an LQR design. 
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The result given in the Theorem 5.1 is extended in the 
Theorem 5.2, to all perturbations in system matrix A and output 
matrix C of the form XAA and XAC where X is a scalar such that 0 < 
X < 1 

The perturbed system equation now can be written as 
x ( t ) = (A + X AA) x(t) + B u(t) (5.1a) 

y(t) = (C + X AC) x ( t ) . (5.1b) 

Remarks 

The Theorem 5.1 is useful when the perturbation matrix AA is 
fixed. The above result have applications in the design of robust 
decentralized PI controller for LT I interconnected system by 
treating interconnection matrix as a perturbation matrix AA. It 
also has applications in the design of robustopt imal PI 
control 1 et s fot system with multiple operating points . 

Theorem 5. 2 

If the condition (5.14) is satisfied for any AA and AC by 
the control law (5.8) then It will also be satisfied for the 
perturbation XAA and XAC by the same control law for all 0 5: X> 1 

Proof 


The following inequality will be used for the proof 


If S > 0 and T > 0 
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then XS + (1 - X)T > 0 for all 0 < X < 1 (5.16) 

where £ and T are symmetric: matrices . 

Satisfaction of condition (5.14) gives 

Q + 2otP - TaA 1 P + P AaI is 0 

Ot ^ O t Ot J 

Also 

Q + 2otP > 0 

Ot 

as 0 5: 0 and P >0 

at 

Using (5.16) we obtain 

+ 2<*e a - [aa* T„ * p„ Al]j + [l - a][q ♦ 2«pJ > o 

= > x£q + 2c*pJ - JxAA* XAaJJ >0 (5 17) 

or 

XAA 4 P^ + KAA - + 2c*pJ so (5.18) 

Proved . 

Remar ks 

The Theorem 5.2 Implies that if the optimality of a closed 
loop system using optimal PI controller is preserved for some 
perturbations AA and AC in state matrix and output matrix, then 
the optimality will also be preserved for perturbations XAA and 
XAC, where 0 ^ X < 1. 
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It is particularly useful when only one element of either A 
or C is subject to perturbations or the perturbations are linear 
function of a parameter. In such cases, results of Theorem 5.1 
and Theorem 5 2 can be combined to find the range of the parameter 
variation which will ptesetve the optimality of the perturbed 
system using optimal PI controller designed for the nominal 
system. 

Now a modified version of the Theorem 5.1 is given which 
gives a condition which when satisfeied guarantees a reduced gain 
margin and reduced phase margin than guaranteed by optimality 
condition for perturbed system using optimal PI controller 
designed for nominal system. 

Theorem 5. 3 

If 

T 

F = r T a - B K 1 P + P Ta -Ik]] + ^ R K (5.19) 

« LI p 04 « l p <*J J « « 

(where ft > 1 , is a scalar) is atleast n.s.d. for some a > 0 and ft 
> 1, then the perturbed system (5.4) using the optimal PI 
controller designed for optimal system, will be robustly stable 
with reduced phase margin of cos * (>- -^), reduced gain reduction 
tolerance of ^ % and with infinite gain margin. 

Proof" 


The proof follows on the same lines as given by Hole and 


Mahanta [11] for an LQR design. 
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I£ the perturbations AA and AC are not fixed and also not of 
the form XAA and XAC, 0 < X < 1, then the following theorem gives 
useful condition, the satlfaction of which preserves the 
optimality of the perturbed system. 


Theorem 5. A 


If the perturbation matrix AA satisfies the condition 


AA 


* P. 


X . (Q + 2«P ) 
mm or 

2X (P ) 
max or 


(5.20) 


Then the pertuibed system will satisfy the optimality condition 
(5.12). 


Proof 


Starting from condition (5.14), we get 


Q + 2otP a - ^AA* 


P + P AAl £ 0 

ot ot I 


or 


AA* P + P AA £ Q + 2o*P 

Ot Ot Ot 


=> 


max 


|aa* 


P + P AAj < X . IQ + 2<xP 1 
ot ot J nun j_ otj 


(5.21) 


or 


AA* P + P AA 

Ot ot 


S + 2 " P J 


(5.22) 


as j | U | j = X [ U ] where U is a symmetric matrix 
s max 


Since 




form (5.22) it follows 


AA t P + P AA < X . |Q + 2«P 1 . (5 23) 

at I cuf I | min ot 

s s * * 1 s 1 s L J 

Since 

= I AA 
s * 

from (5.23), we get 






Uhile designing a robust controller, a designer usually has 
information of the bounds on each element of the perturbation 
matrices AA and AC. The following theorems give conditions which 
can readily be applied to evaluate the optimality robustness of 
the optimal PI controller. 

Theorem 5. 5 


If the pertubation AA satisfies the condition 
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x . 

ram 


[Q + 2otP ] 

OL 


2 X 

max 


(P ) 

a. 


(5.24) 


then the closed loop perturbed system using optimal Pi controller 
designed for the nominal system will satisfy optimality condition 
(5.12). 


Proof 


The proof follows from Theorem 5.4, since 



E 


£ 



Theorem 5. 6 


If denotes some given bounds on all the elements of the 
perturbation matrices AA and AC, i.e , on all the elements of AA 
Then if 


si 


1 

2-/(n+p)n 


X . 
ram 


[Q +2« P ] 

ot 


A. 

max 


(P ) 

a. 


(5.25) 


the closed loop perturbed system will satisfy the optimality 
condition (5.12). 


Proof 

From the definition of Euclidean norm it follows that 
| |AA| | E £ V (n+p)xn 

(since only (n+p)xn elements are subject to perturbation m AA) . 
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Using the above relation and the result of Theorem 5.5, the 
proof of Theorem 5.6 follows directly 

Remark 


Generally all the elements of AA and AC are not subject to 
variation. If i elements, where i si (n+p)xn of AA are subejct 
to variation then the condition (5.25) will reduce to 




A 


< 


1 

2VT 


x. . 

lam 


CD +2ot P ] 

a. 


X 

max 


(P ) 

CA 


(5.26) 


a h in this case 


l|AA|| E * vr « A . 

It may happen that for some given Q and R there is no ot to 
satisfy the condition (5.20). In such cases a modified condition 
can be used, which when satisfied guarantees a phase margin less 
than 60°, a gain margin of gg and gain reduction tolerance of less 
than 50%. The modified condition Is given In Theorem 5 7 


Theorem 5-7 s 


If the perturbation matrix AA satisfies the condition 


1 1 1 1 . ^ f* 


ab 


X . j 0 + 2otP + P— 0- K t R K 1 

_ mxn [ « ( ft J ca otj 


(5.27) 


2 X (P ) 
max or 


then the perturbed system using optimal control law (5.8) will 

-I 1 

have a phase margin of atleast cos (1- ^-r) , a 

SO 

infinity and gain reduction tolerance of -r— %. 


gain margin of 
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|AA|| 


mm 


Q + 2oiP + 

oi 


ab 


M 


K 1 E K 

« ot I 


(5.28) 


2X (P ) 
max or 


then the perturbed system using the optimal control law (5.8), 

-1 1 

will have a phase margin of atleast cos (1 - -^ 5 ) » a gain margin 
of infinity and gain reduction tolerance of % . 

Proof 


Proof follows from Theorem 5.7 on noting that 
1 I*A|| e * ||AA|| s . 

Theorea 5. Q 


I f denotes some given bounds on all the elements of 
A 

perturbation matrices AA and AC, i . e . , on all the elements of AA 
Then if 


ab 


1 

2 Y(n+p )n 


F_ — (73-11 _t — 1 

X . Q + 2«P + K R K 

min { ot ( ft j a 

X (P ) 
max o. 


(5.29) 


the closed loop perturbed system using the optimal control law 
(5.8) will have a phase margin of atleast cos * ( } - i-r), a gain 
margin of ct> and gain reduction tolerance of ^ %. 

Proof" 


From the definition of Euclidean norm it follows that 
||AA|| e £ V(n+p)xn « A 

(since ony (n+p)xn elements are subject to perturbations in AA) . 
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Using the above relation and the result of Theorem 5.8, the 
proof of Theorem 5.9 follows directly. 

Remarks 


If l elements where l < (n+p)xn of AA are subejct to 
variation the condit ion(5 29) will reduce to 


ab 


1 

YT 


mi 


at 


Q + 2«P 


,P ♦ ffcil K l E K 1 

p 1 ° "I 


ot 


(5 30) 


X (P ) 
max cn 


As ft — ► oo, the phase margin tends to zero. Hence the system 
will be at the verge of Instabi 1 i ty , though the gain margin 
remains go. 


5. 4. NUMERICAL EXAMPLE 

This is an example of an aircraft considered by Ulnsor and 
Roy [25]. The longitudinal motion of the aircraft is described by 
the 4th otder system equation. The vailable state matt lx is 


A+AACIl^) 


0 


1 


1.401x10 
-2 . 505x10 


M -.558 

q 

1 


-.561 


0 


0 

0 

-1 . 9513 

. 0133 

-1 3239 

- 0238 

.358 

-.0279 


The nominal value of M Is -1.48 and the input matrix B is 


constant . 


B = 


0 

-5.3307 

-.16 


0 

6 .447x10" 
-1.155x10 


-.2669 
- .2511 


0 


. 106 


.0862 
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The output matrix C la alao assumed constant 


1 . 00x10 


1.13 


1 12 


2x10 


5x10 


The augmented system will be a system of order 7. The augmented 
state matrix and augmented input matrix wil be given by and Ij, 
where 

[0 1 0 OOOOl 


1 401x10 


-2.505x10 


11 - 558 -1 . 9513 0133 0 0 0 

q 

1 -1.3239 -.0238 000 


A =A+AA 
P 


- 561 


1x10 


1 .12 


.3580 -.0279 000 


.23 000 


1.13 


2x10 


5x10 


.11 000 


.1 0 0 0 


-5 3307 


6.447x10 


- .2669 


-.16 


-1 155x10 


- .2511 


. 106 


,0862 


Taking Q = 1^ and R 


I „ , when M takes nominal value, l.s , 
3 ' q ' ’ 


when M = -1.48, AA = [0] When M = -4.492 

q q 




and when M = -.032 

q 

0 0 0 0 0 0 0 ■ 

0 1.448 00000 

0 0 0 0 0 0 0 

AA = 0 0 00000 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

For these two values of AA , i . e . , corresponding to M - -4.4 92 and 

q 

H = -.032, the optimality condition (5.13) is satisfied for ot = 
Q 

0.1 Using the result of Theorem 5.2, we can say that the 
condition (5.13) will be satisfied for all XAA, 0 < A < 1, for ot = 

0.1 When X = 0 , the perturbation m system matrix will be zero 

corresponding to M = -1.48. Uhen X = 1, the perturbation is AA 

corr esponding to M = -4.492 on one extreme and M = - 032 on the 

q q 

other extreme Therefore for all -4.492 < M < 032 the 

q 

perturbed system will satisfy the optimality condition (5 13) for 
ot = 0.1. For increasing this range of for the same degree of 
stability, i.e., for ot = 0.1, we will have to sacrifice in terms 
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of phase margin and gain reduction tolerance by taking ft > 1 . For 

a = 0.1 and ft = 1.2 the perturbed system will satisfy the 

condition (5.19) for all -5.117 ^ M ^ .728. Thus the phase 

Q 

~ i 1 o 

margin reduces to cos (1 - Yft^ = S4.2 3 and gain reduction 

tolerance reduces to % = 41 7 % . 

Now the bound on the perturbation matrix AA is calculated 
which will preserve the optimality of the perturbed system. 

From (5.20), 


^a S 


A (Q + 2otP ) 
max « 

2 a (P ) 
max ot 


(5.20) 


— -3 

For Q = I ^ and « = 0.1 the value of comes out to be 3 .599x10 

If the condition (5.20) is not satisfied by this bound , we can 

Increase this bound by sacrificing the phase margin and the gain 

reduction tolerance by taking ft > 1. This increased bound for ft>l 
la given by (5.27). 


mm 


ab 


Q + 2<*P 


ot 


OL 


JL1 


O i 


1 


2 A (P ) 
max or 


-3 

The value of fj . for ft = 1.2 comes out to be 3.633x10 , which is 

ab 

greater than jj calculated earlier. 

5. 5 CONCLUSION 


The optimality robustness of the optimal PI controller is 
discussed. Some sufficient conditions are derived which when 
satisfied guarantee optimality for the perturbed system, with 
perturbation in state matrix A and output matrix C, using optimal 
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PI controller designed for the nominal system. These conditions 
are in terms of the parameter the state weighting matrix Q and 
input weighting matrix R, thus helping the designer to make proper 
choice of «, $ and R. 

Some numerical examples are solved to illustrate the 
application of the optimality robustness results. 



CHAPTER - 6 


CONCLUSION 

In this thesis the design method of an optimal PI controller 
is discussed and its optimality as veil as stability robustness 
properties are discussed. 

The Chapter 2 discusses the PI controller design method 
The design is carried out by converting the PI controller design 
problem into an LQR design problem. The optimal control lav thus 
obtained is a linear function of states and the integral of error 
between the output and the reference signal 

In Chapter 3 and Chapter 4 the stability robustness of the 
PI controller is discussed It is shovn in Chapter 3 that the 
optimality of the systems using the optimal PI controller 

O 1 

guarantees a phase margins of atleast 60 and a gam margin of 
to to. Some bounds on the perturbations In the system matrices are 

obtained in Chapter 4 £ ot the stability of the perturbed system 

using optimal control law designed for the nominal systems. 

The optimality robustness results are discussed In the 
Chapter 5, where sufficient conditions are derived which when 
satisfied guarantees optimality for the perturbed system. Some 

bounds are also obtained on the perturbations In the state matrix 
which preserves the optimality of the perturbed system. 

Numerical examples ate solved at the end of Chapters 2,3,4 


and 5 to illustrate the results. 
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APPENDIX - 1 

Starting from the ARE (2.9) 

A* P+ PA + Q - P B R _1 B fc P + Q = 0 (A 1.1) 

O O O O \ s 

Add and aubttact aP £t om (At.l) to obtain 

o 

^-sl - aJ P o + P o ^ si - A J + P o B R _1 B t P o = Q (A1.2) 

From equation (2.8b) 

K = R -1 B t P (A1.3) 

o 

Putting R ^ B* = K into equation (A1.2), we get 

f-sl - a] P + P fsl - a] +K t RK = Q (A1.4) 

i J ° °l J 

Multiply on the left by R ^ 2 B* (-si - A t ) * and on the right by 
(si - A}" 1 B R _1/2 to get 

R~ 1/2 B t [-si - A*j P q B R~ 1/2 + R~ 1/2 B t P q [si - aJ B R~ 1/2 
+ R _1/2 B* £-s I - A*J K* R K q [si - aJ B R~ 1/2 

= R _1/2 B* ^-s I - A*J Q [si - aJ B R~ 1/2 (A1.5) 

Add I to both side of (A1.5) 

Since R -1 I* *o ~ * 

=> R -1/2 I t P = R 1/2 K 

o 
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or P o I R" 1 / 2 = R 1 ' 2 

Using above relations and replacing s by j* W6 get 

[’ * Rl/2 - 5 f B «- 1/2 ][: * * 1/2 Kf^x _ X ]- 1 S ,-*/*] 

1 + R B (~ J “ I ' At ] S [ JttI - a] B R~ 1/2 (A1.6) 

o is P-s.d, we get following inequality from (A1.6; 

[ I + Rl/2 ~ *)' * *" 1/2 ] t [l * K - A^B R-/ 2 ] , X 

(A1 . 7) 
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APPEMDIX - 2 


For obtaining the gain margin and phase margin results 
discussed in Chapter 3 of this thesis following r wo lemmas are 
useful . 


Lemma A 5 Let V, U be seqtiare 
dimensions such that 

+ V*J f 1 + vj * I 

U* + U > I 

& 

then I + VU is non singular ([.] 
transpose of matrix [.]). 

Proof* 


omplex matrices of the same 

(A2 1) 

(A2.2) 

•epresents the complex conjugate 


Suppose (I + W) u = 0 for some u & 0 


= > VUu = - u 


Now + V*J (l + V ] ~ 1 


=> I + V* + V + V* V > I 


= > 


* f * 

U V 


+ V + V V 


j B i 0 


■ > U* V* U + U* V U + W* V* V U > 0 


= > u* Ju* V* U + U* V U + U* V* vuj u > 0 
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Putting V W u = - u 

* * * 

and - u = ti U V 

We get 

# A * A 

-u Wu-u U u + u u £ 0 


= > 


u 


' >1 
U 


W - I 


u 


= > 


U + W - I 


=> U* + U < I 

This is a contradit ion , hence I + V U is nonsingular . 

Lemma B s Consider the closed loop system def ined in Fig A2 . 1 . 

Suppose there are no unstable pole-zero cancellations 
in the product V W, and that for all V( j<o) and U( j<o) 
satisfy (A2.1) and (A2.2). Suppose also that with W 
replaced by I the closed loop is stable. Then the 
closed loop of Fig . A2 . 1 is stable . 

Proof" 

Replace U by = (1 - ■£ ) I + <=W where <s can vary in [0, 1]. 

Observe that <e = 0 corresponds to a known stable situation. 


and e 


1 to the situation of interest. Further 
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»* + = 2 (1 - e) I + « (U* + B) 

> 2 (1 - €) I + el 

= (2 - <=) I 

Hence 

U* + U > I 


Also, since (I+V)(I+V)>I, by assumption, then from lemma 

A, I + V(jo>) U ( j<o) is nonsmgular for all co, and for all <= 
€ 

between 0 and 1 . 

Now with U t eplaced by U , the closed loop transfer function 

€ 

matrix is VU (I + VU ) * As -e varies from 0 to 1 , an instability 
€ € 

can arise only when a closed loop pole moves from the open left 
half-plane to the right half plane. In so doing it must cross the 
axis. That is I + V(j^) U ( j*o) becomes singular for some « 
belonging to [0, 1] and some which is a contradiction Thus 

there is closed loop stability with U present. 


The above two Lemmas can be combined to say that when 


[ r + v *) i 1 * v ] 

and the system in Fig. 
in Fig. (A2.2) will be 

U* + W > 


> I 

(A2.1) is stable for U 
stable when 

I . 


I then the system 
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APPENDIX - 3 


Pr oof of the The or em 4. 3 

This has been taken from [17] It is given here for the 
sake of completeness. Consider a system whose dynamics are given 
By 

x(t) = A ffl x(t) + B ffl u(tl + f [ x(t; , u(t; ] (A3.1) 

where x( t ) e Q? n , u( t ) <s tR 10 . 

A and B denote the system model parameters and f is a 
mm r 

nonlinear vector function. 

The nominal system is given by 

x(t) = A ffl x(t) + B ffl uC 1 1 (A3.1) 

For designing LQR the performance index to be minimised is 
<o 

J = j e 2oit [x T (t) Q x(t) + u T (t) R u(t)] dt (A3. 3) 

o 

where Q £ 0 and R > 0 

1/2 

Under the assumption of (A , B ) controllable and [A , Q ] 

fid in m 

observable the optimum control comes out to be 

u(tj = - R -1 B* P x(t), (A3. 4) 


where P > 0 is the solution of the following ARE 
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Ta + oil] P + P [a + oil] - P B R~ 1 B* P + Q = 0 (A3. 5) 

L m J l m J m m 

The resulting closed loop system Is given by 

x(t) = (a - B R -1 B* p] x(t) + f (x(t) ) (A3 . 6 ) 

L in m m j 

Tiheo rem i Let D = Q + PB R ^ B* P , then if the nonlinear 

m m 

vector function £(x(t)) satisfies the condition 

J j f(x(t){ j min X(D) ot min X(P) 

^ ^ = + 

| j x( t ) j | 2 max X(P) max X(P) 

the closed loop system (A3. 6) Is stable. 

Proof* of Lho theorem s 

The Lyapunov function Is chosen as 

V(x) = x T P x (A3. 7) 

where P ia the solution of (A3. 5). 

Ue require V(x) :£ 0 for the stability of the closed loop 
system (A3. 6). Substituting (A3. 6) into 

* * T T * 

Vlx) = x Px+x Px 

and writing A for A -B R ^ B t P, we obtain 
& c mm m 

V(x) = x T Ja* P + P Aj x + 2 f t (x) P x (A3. 8) 

Substltion of (A3. 5) into (A3. 8) yields 
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V(x) = - x^ D x - 2«x^ P x + 2 £* (x) P x 

Since P > 0, 

£*(*) P * s ||f(x)|| E | |Px| | E s ||£<x)|| b | | P| | g | |x| | e 

r mm A.(D) ot mm A-(P ) 1 

* -{ + 1- ®ax k ( P ) j | x j | f 

t 2 max x(P) max x(P) J 

The above inequality is due to the assumption in the theorem. 
=> f t (x) P x ■£ { b “in D) + « min a^(P) i- jjxllp 

(The above inequality holds when the inequality in the theorem is 
satisfied;. Thei e£or e , 

I 1 + 2 at 

fl- 
it Is easy to see that , for ot > 0 , V(x) ^ 0 for all x(t) and 

hence (A3. 6) Is stable. 

Thus the theorem Is proved. 


( p - 


min A. ( P ) 


V(x) = 


‘I-- 


mm A.( D) 



